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Abstract 

We have formulated a theory of self-assembly based on the notion of local 
gauge invariance at the mesoscale. Local gauge invariance at the mesoscale 
generates the required long-range entropic forces responsible for self-assembly 
in binary systems. Our theory was applied to study the onset of mesostructure 
formation above a critical temperature in estane, a diblock copolymer. We 
used diagrammatic methods to transcend the Gaussian approximation and 
obtain a correlation length £ ~ (c — c*)~ 7 , where c* is the minimum concen- 
tration below which self-assembly is impossible, c is the current concentration, 
and 7 was found numerically to be fairly close to 2/3. The renormalized diffu- 
sion constant vanishes as the critical concentration is approached, indicating 
the occurrence of critical slowing down, while the correlation function remains 

finite at the transition point. 
61.25.Em, 64.60.-i 
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I. INTRODUCTION 



Microphase separation is the tendency in certain mixtures such as amphiphilic fluids and 
diblock copolymers, of one of the components to form mesoscale aggregates of the size of ~ 
100 A. Such self-assembly is to be distinguished from the phenomenon of nucleation in single- 
component fluids. Nucleation is the precursor of a phase transition, and as such indicates 
an instability. Self-assembly on the other hand indicates the ability of a given mixture 
to grow islands of one of the components to the size of a couple of hundred Angstroms, 
and then stabilize the growth, so that phase separation occurs only on a mesoscale, rather 
than a macroscale. Above a certain temperature T*, mesoscale structures are formed in a 
random fashion, as long as the concentration of the self-aggregating component is greater 
than some minimum value c*. As the concentration is increased continuously above c*, 
the self-assembling systems first form spherulitic structures, changing to fibrillar and then 
lamellar structures.0 This is inferred experimentally using small-angle (X-ray or neutron) 
scattering. The correlation functions in k-space obtained from such experiments display a 
peak around some wave- vector indicating the average spacing between these islands.!'! The 
width of the peak represents the spread in the average spacing of these islands. Below T*, 
the mesoscale aggregates form regular arrangements (e.g. hep, fee, etc) via a first order 
transition.!!! The regularity of these lattices can be inferred from small-angle scattering 
experiments which display harmonics of the main peak.! 

It is generally believed that self-assembly in mixtures is due to the competition between 
the tendency of the components to phase separate on a macroscale, and a long-range entropic 
(statistical) force caused by the presence of chemical bonds linking the components in the 
mixture!] In the case of amphiphilic mixtures, it is the surfactant molecules which provide 
the glue which allows mesoscale segregation to occur. In the case of diblock copolymers, end- 
groups on the two species create inter-specie bonds, thereby playing the role of a surfactant. 
A molecular-level description of mesoscale structures (micelles) in liquids, and aggregates in 
copolymers that are a couple of hundred angstroms in size is a challenging problem. At this 
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scale, raw simulations which begin at the molecular level are simply impossible to perform 
for realistic molecules with the current computational technology. 

Parallel field theoretic efforts in both amphiphilic fluids as well as diblock copoly- 
mers have been developed over the years to provide an understanding of microphase 
separation.!IH3 We will show in this paper that the principle of local gauge invariance 
with respect to the £0(2) group can be applied successfully to unify the above theories 
with a common thread, and furthermore, to derive a generalization of these mesoscopic 
theories of self-assembly.0 We have interpreted the gauge fields we obtain as giving rise 
to statistical correlations between concentration fluctuations. These statistical correlations 
could be thought of as effective interactions which arise at the mesoscale from the underly- 
ing Coulombic interactions at the molecular level, between the components of the mixture. 
While the use of local gauge invarianceEl is quite well established in particle physics, its 
usefulness in settings other than quantum field theory (QFT) is appreciated only under rare 
circumstances.^ We note that while the dynamical use of local gauge invariance is novel at 
the mesoscale, gauge theory has been used routinely in the past to classify defects in con- 
densed matter physics.0 Our theory is applicable to diblock copolymers, oil-water-surfactant 
mixtures, and in general any self-assembling system, e.g. binary alloysJi! 

A further importance of our paper lies in the fact that we have gone beyond the Gaussian 
approximation or the Mean Field Approximation (MFA) used conventionally in meso-scale 
theoretical investigations.!'! While the MFA may be a reasonable approximation to study 
self-assembling systems far from phase transitions, it is obvious that one must necessarily 
go beyond the MFA or the Gaussian approximation in order to properly study the on- 
set of self-assembly. To be more precise, we point out that investigations of the onset of 
self-assembly (as the composition is varied) in the literature&iS yield a correlation length 
diverging with the square-root signature of the MFA. Experimental observations cited by 
Woo et a 10 suggest that the true exponent is larger than |. We will take seriously in this 
paper the quantitative suggestion of Woo et al that there is a need to go beyond the MFA 
or the Gaussian approximation to study the onset of self-assembly. 



There have indeed been investigations in the past where Renormalization Group (RG) 
techniques have been used to study the first order transition from a disordered to an ordered 
phase in copolymers, as the temperature is varied.@@~00 But the onset of the self-assembly 
of mesoscopic structures into a random arrangement above a critical temperature, as the 
concentration of one of the components of the binary mixture is varied, is an issue that has 
not been addressed theoretically in much detail beyond the Gaussian approximation. Our 
investigation reveals that this transition is analogous to the critical point in phase transition 
theory, in that the correlation length diverges as a 2/3 power law, and the diffusion constant 
goes to zero, implying critical slowing down. But the correlation function itself does not 
diverge. In this sense, we are investigating a Lifshitz point .111 Our detailed calculations apply 
specifically to estane, a diblock copolymer. However, one may invoke universality arguments 
to argue that our results are applicable more generally. 

Our theory represents a generalization of the 4 field theory proposed by Landau and 
Ginzburg to study phase transitions. While our theory is slightly similar to the standard 
Landau-Ginzburg theory in that they are both nonlinear and deal with an order parameter, 
it is clear that there are some major differences. First of all, our theory is nonlocal in 
character. Secondly, the nonlinear term in our theory not only contains a cubic term (in 
addition to a quartic term), which arises naturally from an expansion around the average 
value of the fields, but the nonlinear term also contains derivatives of the concentration. 
The derivative form of the nonlinear coupling is dictated by the fact that ours is a gauge 
theory, in which covariant derivatives are defined. We have used this theory to investigate 
the onset of self-assembly in estane, a diblock copolymer, and we found that the correlation 
length diverges with a power which is fairly close to the universal value of (2/3). We also 
found that the renormalized diffusion constant goes to zero as the minimum concentration 
c* (below which self-assembly is imposible) is approached. 

We foresee a rich variety of applications of our approach to other questions regarding 
diblock copolymers, such as their viscoelastic properties. We also foresee investigations of 
time-dependent phenomena in copolymers, such as detailed studies of critical slowing down 

4 



at the onset of self-assembly. 



II. THE GAUGE THEORY 

The starting point of our mesoscale theory is an internal energy functional which is 
quadratic in the gradient of a two-dimensional vector. For the moment, we will consider 
isolated systems, so that the quantity that is conserved is the internal energy.il We will 
shortly consider entropy effects as well. Consider the following form for the energy functional: 

PU = pJuo(c(s))d 3 s (1) 

n t , w f 9\ <9c*(s) <9c(s) , . 

/*.(«(.)) = (f) -^-^ (3) 

where t indicates a transpose, repeated indices are summed over, and, 

In the above equations, s is a dimensionless co-ordinate variable, k is Boltzmann's con- 
stant, T is the temperature, Ch is the number concentration of the first type of specie, and 
c s is the number concentration of the other specie in a binary mixture. The concentrations 
are normalized to the total number concentration. The constant g is essentially a dimen- 
sionless diffusion constant. Such energy functionals have been considered over many years 
as contributing to the total internal energy of binary mixtures. i'0 We will use this from as 
our starting point to generate a more complete energy functional using gauge invariance. 

From Eqns.(l)-(4) we see that u Q is invariant under global rotations of the vector c. 
These are rotations in two dimensions, and the appropriate group to consider is 50(2). The 
physical origin of this group can be traced back to the fact that the quadratic (positive, semi- 
definite) form of the energy density (Eqn.(3)) is dictated by expanding the internal energy 
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around a minimum, in a Landau-like fashion. The form of the energy density contains 
gradient operators, which permits us to perform 5*0(2) transformations around not just the 
origin in (ch, c s ) space, but around any arbitrary fixed vector in this space. In particular, 
we shall use 5*0(2) around the vector defined by the average concentration of each specie 
viz., (c° h) c°). This is a natural representation for our system, since our final goal is to study 
self-assembly in binary systems, characterized by local, mesoscale fluctuations around the 
average concentrations. 50(2) transformations of these fluctuations demands that d\ + c' 2 s 
= constant, where c'h and c' s denote deviations of the specie concentrations around their 
averages. Thus 50(2) transformations can cause the components of (c'h,c' s ) to become 
negative. But this is acceptable, since concentration fluctuations around the average can 
indeed be negative or postive, as long as the total concentration for each specie does not 
become negative (see Eqn.(ll)). In what follows we shall be tacitly performing local 50(2) 
transformations around the average concentration vector (c£,c°), culminating in Eqn.(13), 
which is a central result in our paper. 

Our physical motivation for seeking local gauge invariance of d\ + c' 2 s under 50(2) is 
the same as that of Yang and Mills0, and in quantum electrodynamics, where one observes 
the invariance of the noninteracting Lagrangian, which is bilinear combinations of the fields, 
under certain global transformations. One then demands covariance of the theory when these 
symmetry operations are local i.e., when the transformations are space-time dependent. A 
reason for this, as given by Yang and Mills, is that one can now freely interchange between 
the fields as one moves through space and time, while leaving the physics covariant. It is 
intuitively clear that such is the case in our problem, where chemical connections between 
the two species n our system allows for an admixture of the two components, rather than 
permitting a complete phase separation to occur on a macroscale. Thus, instituting local 
gauge invariance under 50(2) in our binary mixture is equivalent to allowing interactions 
between the components. Beyond this initial motivation, it is equally important to show 
that the result of local gauge transformations of u a lead to physically significant results as 
epitomized by Eqn.(13). 
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We remark in passing that u Q is also invariant under the translation group T2, where we 
consider the transformations c — > c + a. Based on the work of Edelenlll in solid mechanics, 
we believe that seeking local gauge invariance of u Q under T(2) may lead to a study of defects 
in our system. 

Following Yang and Mills,0 local gauge invariance of u Q under 5*0(2) motivates us to 
define new fields b, which have invariance properties appropriate to SO(2). We define a 
covariant derivative — > (^- + qrbi), where r is the generator of SO(2), q is a 'charge', 
or equivalently, a coupling constant, and the 6-fields are analogs of the magnetic vector 
potential in electrodynamics. These 6-fields give rise to effective interactions between the 
hard and soft segments of estane. These effective interactions are to be thought of as arising 
form the underlying electrostatic interactions between molecules, monomers, etc. The energy 
functional for the 6-fields is defined a la Yang and Mills, via the minimal prescription. With 
this, our original internal energy density is transformed into: 

(3u D — ► f3u = (3u + (3u int + (3u Y m (5) 

where u int refers to the interaction energy density, and uym is the energy density as- 
sociated with the Yang-Mills 6-fields alone. Equivalently, we may define the total energy 
functionals associated with these energy densities: 

PU ^(3U = (3U + (3U mt + PU YM , 

where 

(3u int = Ji(c)bi(s) + 6i(s)/(c)6i(s) (6) 

with 
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r is given by:E3 

-1 

1 

From the above equation, it can be shown that 



(9) 



x dsi 

We need one more definition for completeness: 



Ji(c) = qg ( %^rc(s)) (9a) 



,1\ (dh dbA (dh dbA 



This equation can be cast into the following form: 

Pu Y m = -(j^ k^X (10a) 

Eqn.(lOa) is obtained via an integration by parts, in the transverse gauge. Since we 
are dealing with an Abelian gauge theory, it is permissible to insert this transverse gauge 
manually, without resorting to the formal machinery of Faddeev and Popov. 

Note that we are utilizing a non-relativistic version of the Yang-Mills procedure, since we 
are only concerned with time-independent problems. Furthermore, since we are concerned 
with rotations in two-dimensional space, there is only a single generator for the group 50(2) 
(see Eqn.(9)), so that the resulting functional is only quadratic and not quartic in the ^-fields. 

It is important to emphasize that the usual application of the Yang-Mills procedure in 
QFT implies the existence of fundamental interactions. In our case, we are applying the 
principle of local gauge invariance at the mesoscale. Consequently, we do not expect to 
discover any new fundamental interactions by using gauge invariance. Rather, we interpret 
the new 6-fields as yielding correlations between the concentration fields. As proof of this, 
we will show shortly that our approach leads to a generalization of the theories of Stillinger 
and Leibler, where correlations were invoked on physical grounds to describe mesoscale 
structures. These correlations could also be thought of as effective interactions, which arise 
at the mesoscale from the underlying electrostatic interactions between molecules. We will 
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not address the question of how one can make a connection with molecular scale properties 
in this paper. 

Our approach is analogous to the Landau-Ginzburg theory of superconductors in mag- 
netic fields.il In the Landau-Ginzburg theory, the energy functional involving a complex or- 
der parameter is gauged with respect to the £7(1) group. This permits a successful treatment 
of a superconductor in a magnetic field, and even permits a classification of superconductors. 
We have gone further in our theory, and invoked gauge invariance to study correlations that 
develop at the mesoscale. In our theory, there is no external magnetic field to consider. 

The partition function we need to evaluate is now: 

Q = J \\Vc a 9 (c a ) I] Vb k exp - (3{U + U mt + U YM ) (11) 

a=h,s fc=l,3 

Equation (11) is a functional integral, where the step functions denoted by 9 imply that 
we must restrict integration to positive semi-definite values of the fields. 

Since the 6-fields appear only quadratically in the above functional, it is straightforward 
to integrate over them, and obtain an effective internal energy functional involving only c. 
The result is:0 

pU eff = (3U + (3AU eff = PU -\J dh J dV J,(c(s)) ^ j—L-^ J J,(c(s')) 

(12) 

Note that in doing so, we have ignored an overall trivial normalization constant that 
appears in the evaluation of the partition function Q. This is permissible, as this factor 
cancels during the evaluation of averages of observable quantities. 

To see the connection between this rather complicated functional and the older theories, 
we expand the second term on the right hand side of Eqn.(12) around the average concen- 
trations of the two species (c° and c°) that appear in our theory, and retain only quadratic 
terms. The result is: 

f3U eff ~ (3U - Q J d*s c'*(s)5c'(s) + Q / d 3 , / dV c%)S eXp{ -^- S ' l) c'(s') 

(13) 
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where f2 = g 2 q 2 , g = (|J g 2 ((c°) 2 + (c°) 2 ), T = q 2 g 2 g, and the primes on c denotes 



deviations of the specie concentrations from their averages. The matrix S is defined thusly: 



First of all, we see that by expanding around the average value of the fields, we are in 
essence considering the effect of local SO(2) invariance on the correlations which develop 
between the fluctuations of the fields. Secondly, we notice that if we set g — > 0, we recover a 
model very similar to that of Stillinger and Leibler. Q is the Flory-Huggins parameter, and 
represents the immiscibility of the two components in our mixture. It is prescribed auto- 
matically via gauge theory, as long as g is known. The non-local term in Eqn.(13) gives rise 
to correlations which tend to counteract the effect of Q. This frustration is responsible for 
the formation of mesostructures. From the definition of S, we see that we have retained in 
our model the notion of (pseudo) electroneutrality emphasized by Stillinger! and Chandler et 
a/.i Finally, we note that in general g is not zero, so that we have a screened Coulombic cor- 
relation appearing in the second term of Eqn.(13). In this sense, Eqn.(13) may be viewed as 
being similar to the random phase approximation (RPA) applied to the full functional given 
by Eqn.(12). Note that for small deviations (c' h ,c' s ) from the corresponding concentration 
averages, the step functions of Eqn.(ll) have a negligible effect. 

Equation (13) is one of the main results of our paper. It shows that Leibler and Stillinger's 
theories may be understood in the context of gauge theories. Equation (13) gives credence 
to the notion that gauge theoretic ideas may be valid at the mesoscopic level. 

Equation (13) was obtained by expanding fields around their average values. In this 
sense we have broken the symmetry of our system. Combined with gauge invariance, then, 
we get a Yukawa-type screened potential. This effect may be interpreted by saying that the 
gauge fields have acquired a mass. In this sense, Eqn.(13) exhibits the Higgs phenomenon.il 

Before we can compare our theory with experimental data, we need to consider the fact 
that our system is not really isolated, and may be in contact with an energy reservoir, 
perhaps as it is being acted on by mechanical forces in a stress experiment. For a system 




(14) 
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in contact with an energy reservoir, the quantity that is conserved is the Helmholtz free 
energy!!] A = U — ST, where S is the entropy of the system. The entropy of our system will 
be written in the usual form: 



with this, our theory is formally complete. 

While our theory has been able to reproduce the older theories of self-assembly of Leibler, 
Stillinger and Chandler, we believe that the importance of our approach lies in the fact tht 
it provides a natural way to go beyond the Gaussian approximation and the MFA used 
conventionally in meso-scale investigations. By this we mean that our effective functional 
can be expanded in an infinite series beyond the Gaussian approximation. While the MFA 
is a reasonable approximation to study self-assembling systems far from phase transitions, it 
appears obvious that one must necessarily go beyond the Gaussian approximation in order 
to properly study phase transitions, e,g., the onset of self-assembly. 



To see what lies beyond the Gaussian approximation, it is convenient to invoke incom- 
pressibility, so that we can cast the Helmholtz free energy solely in terms of the concentration 
of Ch the concentration of one of the species in our binary system (with the average concen- 
tration of that specie subtracted from it). It is important to point out that the condition 
of incompressibility is to be imposed after the starting functional Uq has been gauged. The 
condition of incompressibility is to be accounted for during the evaluation of the partition 
function for the system. 




(15) 



III. GOING BEYOND THE GAUSSIAN APPROXIMATION 




(16) 



/3a(s) = /3a (s) + /3Aa(s) 



(17) 



(3 A = (3A + (3 A A 



(18) 
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U D = J d 3 s g(Vc h (s)) 2 (19) 
/3AU eff = - (ipj J dV| d 3 S V^(s') • (g ( 1Ch + c 2 h ) g ) s/s V s c h (s) (20) 
-~ = J d 3 s (c,(s)ln(c fe (s)) + (1 - c,(s)ln(l - c h (s))) (21) 

where (3 = 1/ksT, and 

(l/(c»(l - <■»))) c»(s) 2 + S (Vc(s)) 2 
-QiW + (^)J ^-^K w (22) 

w a J* d 3 sV s c h (s) ■ (jc h (s) + c h (s) 2 ) V s c h (s) (23) 

where Q = g 2 q 2 , g = (j ) g 2 ((c°) 2 + (c°) 2 ), and T = q 2 g 2 g, g is essentially a dimensionless 
diffusion constant, and q is a pseudo-charge that arises out of our gauge theoretic consider- 
ations, c° and c° are the average concentrations of the two individual species in our system, 
and 

9o=(f(c h (s))- 1 -V 2 y 1 (24) 
/Ms)) = \gq 2 {c 2 h {s) + (1 - c h (s)) 2 ) (25) 

a = 4((c°) 2 + (c°) 2 ) (26) 
7 = 4 (c° - 1/2) (27) 

The competition between the Flory-Huggins separation parameter Q and the attractive 
non-local term in a (Eqn. (|22|) ) gives rise to the formation of mesostructures. In obtaining 
Eqn. (|23|) , we have ignored terms linear in the fields, and constant terms, as they do not 
contribute to the density- density correlation function. In Eqn.^3|, we have looked for higher 
order corrections to AU e ff. We have ignored cubic and quartic contributions which come 
from the entropy term, as our diagrammatic estimates indicate that they are negligible. 



The form of AA in Equation 23 is a local form. The local form is obtained by retaining 



only the lowest order non-linear terms in an expansion of the full nonlocal form of the 
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interaction term. It is a reflection of the fact that the full form of the nonlocal interaction 



term is screened on a length scale l/y/2g. The asterisk on the integral in Equation ^3 
indicates that a cutoff in momentum space is to be used in the short wavelength limit, 
kmax = V%9- There would be no need for a cutoff if the full form of AA were to be used. 
We note that Equations and [23| represent a generalization of the 4 field theory proposed 
by Landau and Ginzburg to study phase transitions. 

We will now apply this gauge theory to the onset of self-assembly in estane, a di- 
block copolymer. Estane is composed of hard segments of polyurethane, and soft strands 
of polyester. The hard segments display microphase separation on the scale of (9(100) 
Angstroms. It is appropriate to consider an approximation to the form of a , which ensures 
that Porod's law is satisfied in the small wavelength limit .0 



f3a Q (k) « c£(k) (a + g'k 2 + T'k 4 ) c h (k) (28) 

where: g> = g/(l - l/2((c°) 2 + (c°) 2 ), V = l/(2g 2 ((c°) 2 + (c°) 2 )), and a = l/(2c£(l-cg)). 

The term in parantheses in equation |28| represents the inverse of the structure factor 
(Fourier transform of the density-density correlation function) in the Gaussian approxi- 
mation. Following the formulation of our gauge theory, g' > 0. As such, the Gaussian 
approximation has to be improved upon before seeking agreement with experiment. Unless 
g' gets remormalized to a negative value, the structure function will not yield a peak at 
some non-zero value of the wave number, which would characterize a micro-phase separated 
system. In what follows, we shall drop the subscript h which appears on the field c. The 
partition function is defined as (J is an auxiliary field): 

Q[J] = J Dcexp -p (a (s) + Aa(s) + 2J(s)c(s)) (29) 

We can now use standard perturbation techniques^ to develop a series expansion for 
corrections to the Gaussian approximation. As such, it is meaningful to make certain that 



the dimensionless coupling constant a defined in Equation B7I is less than one. On the other 



hand, it is well-known0 that such series are asymptotic in nature. Thus a < 1 is not a 
panacea. We define the two-point correlation function as follows: 
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.V(x x .x 2 ) = ( - 



* * \nQ[J] 



(30) 

j=o 



5 J( Xl ) 5J(x 2 ) 

Figure 1 gives a pictorial representation of the two terms in the definition of Aa. It 
is clear that each of these two interaction terms yields a separate perturbation series. In 
addition, there will be a series formed out of the cross-terms as well. There are no cross-terms 
up to the 2- loop level. The cubic interaction term first yields non- vanishing contributions 
in second order perturbation theory (Figures 2a-2b). The quartic term yields non- zero 
contributions at the one-loop level (Figures 3a-3b). We have verified explicitly that all other 
(asymmetric) diagrams arising up to the 2-loop level add up to yield a null contribution. 
Similar cancellations are also obtained in theories of dendritic growth.il In our calculation, 
the two series arising out of each of the two interaction terms were evaluated only to the 
first non- vanishing order. 

Figure 2a (tadpole) renormalizes g'\ 

Sg'h) = - ( ' ) a 2 j 2 / dk 7 -— (31) 

where k max = <fZg. 

This tadpole diagram is crucial in helping us achieve agreement with experiment. It 
dominates the contributions from Eqns.^ and^. Note that it has a sign opposite that 
given by Eqn.|33] below. Without EqnfJTl, the renormalized structure function would not 
yield the characteristic peak in scattering data. This tadpole diagram is reminiscent of the 
standard (Hartree) tadpole diagram in many-body physics. The physical importance of this 
diagram is as follows. The bare diffusion constant g, if left unregulated, would tend to 
smooth out concentration gradients in our system. It is the role of the screened " Coulomb" 
interactions, having a statistical origin, which are responsible for self-assembly. And it is 
up to these interactions to counteract the smoothing tendency of the diffusion term. This 
is accomplished as described above, by the tadpole diagram, which renormalizes the bare 
diffusion constant so that the renormalized diffusion constant is less than or equal to zero. 

Figure 2b yields two terms in leading order, one which renormalizes g', and the other 
which renormalizes T': 
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/ 4 \ on t k rnax 



(a + c/'P + PA; 4 )^ 

5Y\l) = - 4)« 7 2 / dk j (32) 

l W io (a + g'k 2 + T'A; 4 ) V ; 



The renormalization of g caused by this diagram is dominated by the contribution of the 
tadpole diagram discussed above. 
Figure 3a (1-loop) renormalizes g': 

/ 9 \ rkmax b 2 

Sg'(3) = ( — )a dk- f- — (33) 

V ; W 2 / io (a + g'k 2 + T'k 4 ) V ; 

Figure 3b (1-loop) renormalizes the constant a in Equation A closed form expression 
for the contribution from this diagram is: 

5a(l) = [-= )a dk- — — — — (34) 

We evaluated the integrals appearing above numerically, and then fitted experimental 
data! on estane. We had three free parameters to manipulate. The three parameters at 
our disposal are g, a dimensionless diffusion constant, q the pseudo charge of our gauge 
theory, and finally the length scale A which we used to turn the spatial lengths in our theory 
dimensionless. It is important to point out that the form of the renormalized structure 
function was crucial in obtaining reasonable agreement with experiment. The results of 
fitting Bonart's data are displayed in Fig. 4. The values needed are listed in the figure 
caption. It is worth noting that the value of the length scale A ~ 100 A we obtained can 
be interpreted physically as the mean distance over which averaging has been performed 
to go from an atomistic description to a mesoscale model. The volume described by this 
length scale can accommodate roughly 10 3 polyurethane (hard segment) monomers. With 
the current values of the parameters, our coupling constant a was just under 0.5. 

Armed with values for our parameters appropriate for estane, we varied the concentration 
c° of the hard segments of polyurethane in estane, decreasing it from Bonart's value of 0.25. 
In this way, we could probe how the location of the peak in the structure factor S(k) 
changed, as we varied the concentration of the hard segments. The location of the peak 
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is a measure of the inverse of the correlation length in the system. The purpose of this 
exercise was to investigate how the correlation length behaves as the onset of self-assembly 
is approached. Figure 5 shows a plot of the inverse correlation length as the concentration 
is varied in the vicinity of c* (the minimum concentration below which self-assembly is 
impossible). We found c* to be approximately 0.1922 while the critical exponent was found 
to be fairly close to 2/3 (0.6542). We also found that as c* is approached, the renormalized 
effective diffusion constant g' R = (g' + 5g'(l) + Sg'(2) + Sg'(3)) — > 0, implying that we have a 
precursor of critical slowing down. Our calculation of this critical exponent is the first which 
goes beyond the Gaussian approximation,!! and makes a prediction for diblock copolymers 
regarding the the correlation length near the onset of self-assembly, and shows that the onset 
is accompanied by critical slowing down. Experiments by Koberstein et alS on polyurethane- 
polypropylene systems indicate qualitative agreement with our theory. The difficulty in 
obtaining quantitative agreement lies in the fact that experimental measurements of the 
structure factor start to be swamped by the direct portion of the beam, and the error bars 
increase dramatically, as the peak marches towards the long wavelength limit. 



IV. CONCLUSIONS 

We have shown that the concept of gauge invariance can be utilized successfully at the 
mesoscale to generate the entropic/statistical long-range forces responsible for self-assembly. 
This was done by using local gauge invariance under the 5*0(2) group to derive the older 
(Gaussian) theories of self-assembly. Our approach allows us to go naturally beyond the 
Gaussian approximation. We computed the first nonvanishing contributions beyond the 
Gaussian approximation to the density- density correlation function from the cubic and quar- 
tic terms in our energy functional. We applied our theory to estane, a diblock copolymer, 
above its critical temperature. We found that as the concentration approaches c* from 
above, the correlation length diverges with a (2/3) power law, and the renormalized dif- 
fusion constant tends to zero, implying critical slowing down. The correlation function 
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however, remains finite at the transition. The divergence of the correlation length can be 
interpreted to mean that as c* is approached from above, the probability of finding mesoscale 
aggregates vanishes, so that the average distance between aggregates diverges. Experiments 
by Koberstein et a@ on polyurethane- polypropylene systems indicate qualitative agree- 
ment with our theory. Universality arguments may be used to argue that our results are 
applicable more generally. 

We foresee a rich variety of applications of our approach to other questions regarding 
diblock copolymers, such as their viscoelastic properties. We also foresee investigations to 
time-dependent phenomena in copolymers, such as detailed studies of critical slowing down 
at the onset of self-assembly, and constitutive relations at non-zero strain rates. 
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FIGURES 

FIG. 1. (a) is a pictorial representation of the cubic term in Aa. Each leg corresponds to a 
factor of c, the field. A di indicates that a derivative of the field is to be taken in the i-th direction. 
A sum over i is understood. The dark circle symbolizes a factor of —cry, the coupling constant. The 
negative sign comes from the argument of the Boltzmann factor, (b) is a pictorial representation 
of the quartic term in Aa. A factor of —a is to be inserted at the intersection. 

FIG. 2. (a) represents the tadpole diagram which is crucial in our calculations, (b) represents 
the setting sun diagram. Both (a) and (b) are second order contributions to the correlation function 
coming from the cubic interaction term, the first order corrections being null. 

FIG. 3. (a) and (b) represent 1-loop contributions from the quartic interaction term to the 
correlation term. 

FIG. 4. Comparison of experimental data and theory for estane. The dark circles are experi- 
mental data points without the background subtracted from them; the solid curve is theory. The 
experimental background is negligible, except for the data point closest to the origin. Both data 
and theory are normalized to the peak value. The parameters used to obtain the fit were: c°=0.25, 
g = 8000^4 2 , q = 0.6x 10~ 3 A~ 2 , and A = 95A The peak in the experimental data and theory 
occurs at about 200 A, indicating the average distance between aggregates. Note that the slight 
peak in the theoretical curve is obscured by the experimental data. It nevertheless exists and can 
be verified numerically. The apparent disagreement between data and theory for large momenta is 
due to the fact that the experimental background, which has not been subtracted from the data, 
becomes important in this regime. 

FIG. 5. This is a plot of the location of the maximum in the theoretical structure factor S(k), as 
the concentration of the polyurethane c° is varied in the neighborhood of the critical concentration 
c*, which is m 0.19216. The open circles denote numerical results obtained in the manner discussed 
in the text. The solid line is the curve given by 0.71406(c° — c *) - 65424 . 
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